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We study the scaling of photon-photon correlations mediated by resonant interactions of pho- 
tons with two-level emitters in a one-dimensional photonic waveguide. Recently a new theoretical 
approach based on the Bethe-ansatz technique has been developed to study transport in an open 
quantum impurity. Here we generalize the approach to study multiple emitters-for example, atoms 
or quantum dots. We derive the exact solution of single and two-photon scattering states, and the 
corresponding photon transmission through the emitter ensemble. Correlations of photons at the 
both ends of the waveguide are examined carefully for one, two and three emitters. We show how 
various two-photon nonlinear effects, such as spatial attraction and repulsion between photons as 
well as background fluorescence can be tuned by changing the number of emitters and the coupling 
between emitters (controlled by the separation). Finally we propose a simple scheme for nonrecipro- 
cal optical transmission in the waveguide by placing emitters with different transition energies. Our 
fully quantum- mechanical approach provides a better understanding of cascaded optical nonlinearity 
at the microscopic level. 



I. INTRODUCTION 

Strong coherent photon-photon correlations at the 
level of few atoms and photons have been demonstrated 
in various quantum optics set-ups [1-14]. The efficiency 
of a single or few atoms to induce strong interactions 
between propagating photons is indispensable to realize 
various logic gates for quantum information processing, 
quantum computation and alternative technologies based 
on switching and amplification functionalities. One inter- 
esting recent proposal to achieve strong coherent photon- 
photon interactions is by confining photons in the re- 
duced dimensions, such as, in a one dimensional (ID) op- 
tical waveguide, and coupling these photons with individ- 
ual emitters in the waveguide [15-24]. Tight confinement 
of light fields in the waveguide directs the majority of 
the spontaneously emitted light from the emitter into the 
transport channels, while local interactions at the emitter 
induce strong photon-photon correlations by preventing 
multiple occupancy of photons at the emitter. Various 
nanoscale systems, such as photonic crystal waveguides 
[25], surface plasmon modes of metallic nanowires [26], 
microwave transmission lines [11, 27], optical nanofibcrs 
[6, 28], semiconductor or diamond nanowires [29, 30] 
would act as ID continuum for photons. Different two- or 
multi- level atoms [1, 4-6, 10, 28], molecules [9], quantum 
dots [26, 29], superconducting qubits [5, 11, 27], nitrogen- 
vacancy centres in diamond [30] are used as an emitter 
to couple with the ID continuum of photons. 

Recently a nonperturbative fully quantum mechanical 
approach based on the Bcthc-ansatz technique has been 
developed to study scattering of photons from an emit- 
ter in ID [15, 18, 19]. The main feature of the approach 
is a proper inclusion of essential open boundary condi- 
tions. This is in contrast to the earlier Bethe-ansatz 
techniques assuming a periodic boundary condition to 
obtain thermodynamic information of the quantum im- 



purity problems. Several interesting phenomena, such as 
a two-photon bound state [15], the background fluores- 
cence [15], a single-photon transistor [16], a few-photon 
diode [18], the scaling of clectromagnetically induced 
transparency with number of photons [19], the stimu- 
lated emission [24] have been studied in different systems 
with a single two-level or three-level emitter coupled to 
a ID continuum of photons using the Bethe-ansatz ap- 
proach. Many of these proposals are already tested in 
experiments [9, 12, 26]. The Hamiltonian describing the 
waveguide-emitter system also elucidates the interactions 
of a light beam with an emitter in three dimension (3D), 
when the beam is designed to mode-match the emitter's 
radiation pattern [8, 9, 24, 31, 32]. Therefore, studies us- 
ing the open Bethe-ansatz method including our results 
in the present paper should be also relevant to a wide 
class of recent experiments [9, 12] in 3D. 

At this point it is interesting to ask how these phe- 
nomena scale with increasing number of emitters in the 
waveguide. It has been demonstrated that non-classical 
photon-photon correlations, such as antibunching and 
sub-Poissonian photon statistics are mostly washed out 
when more than one atom is excited by low intensity, 
monochromatic laser light [33, 34]. The quantum statis- 
tical characteristics of radiation produced by a cooper- 
ative system of two and more emitters in 3D have been 
studied in detail without [35-38] and with [39-41] includ- 
ing the interactions between emitters. Though we are not 
aware of any previous study about quantum features of 
photon-photon correlations for multiple two-level emit- 
ters in a ID waveguide. Here we generalize the open 
Bethe-ansatz approach to multiple emitters. We calcu- 
late the exact single and two-photon scattering states, 
and the photon transmission through the emitter ensem- 
ble. Correlations of photons at the entrance and the exit 
of the waveguide are examined carefully for one, two and 
three emitters. We show how different two-photon non- 
linear effects, such as spatial attraction and repulsion be- 
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FIG. 1. A schematic of a chain of two-level emitters embedded 
in a one-dimensional wave-guide, and the propagating free 
photon modes at the left and the right sides of the emitters. 



tween photons as well as background fluorescence can be 
tuned by changing the number of emitters and the cou- 
pling between emitters. The coupling between emitters 
in the present set-up is determined by the separation be- 
tween them. Later we show that an emitter ensemble 
in a ID waveguide is a suitable model for realizing a 
fully quantum optical-diode at low light intensity. It is 
within the reach of present quantum-optics technologies 
to control the number of emitters in a photonic waveg- 
uide with a strong emitter-photon coupling. Our fully 
quantum-mechanical approach for transport in open sys- 
tems provides a better understanding of cascaded optical 
nonlinearity in ID at the microscopic level. 

In Sec. II wc outline the open Bethe-ansatz approach in 
detail for two emitters in a ID waveguide and propagat- 
ing photon modes. Single and two-photon dynamics are 
investigated respectively in Subsec.IIA and II B. Nonre- 
ciprocal photon transmission for emitters with different 
transition energies is demonstrated in Subsec.IIC. Re- 
sults of single and two-photon dynamics for three emit- 
ters in the waveguide are obtained in two subsections of 
Sec. III. We conclude with a discussion of our main results 
and some future directions for further study in Sec. IV. 



approximation for emitter-photon coupling is given by, 

Ha = -i J dx [vgd\(x)d x ai(x) + Vgdl(x)d x a 2 (x)}, 

H A = fii|2) n (2| + 2 |2) 22 (2| + J{a x _a 2+ + a 2 _a 1+ ), 

H c = (tW&i(0) + tWa 2 (0) + H.c), (2.1) 

where the first term Ho represents the propagating pho- 
ton modes with group velocity v 3 g at the left (j = 1) 
and the right (j — 2) side of the chain; the second term 
corresponds to the chain of emitters; and the last term 
describes an interaction between the emitters and the 
free photons. Here Clj is a transition energy of the jth 
emitter and J is a coupling between the emitters. V\ (V 2 ) 
is a coupling between the photons at the left (right) side 
of the chain and the emitters. a\(x) and a\(x) denote 
creation operators of a photon at the left and the right 
of the chain respectively. <jj- = c\ g Cj e (<Jj + = cJ- e Cj S ) 
is a lowering (raising) ladder operator of the jth emit- 
ter where ct (ct e ) is a creation operator of the ground 
state \l)j (excited state |2)j) of the jth emitter. Also 
one can include spontaneous emission from the excited 
states of the emitters to other photon modes, which are 
not part of the propagating free photon modes in the 
above Hamiltonian, by including an imaginary part to 
the energies of the excited states of the emitters within 
the quantum-jump approach [19, 43, 44]. Next we scale 
the free photon operators to absorb the group velocity, 



and redefine yv g a,j(x) = cij(x), Vj/yVg = Vj. There- 
fore we rewrite Hq = —i J dx^- =1 2 a,j(x)d x aj(x) and 

He — 5Zj=i 2 Vj a j+ a j(®) + H.c. Hereafter we always 
consider the transformed Hamiltonian. 



II. TWO EMITTERS IN A ID WAVEGUIDE 

We consider a chain of two-level emitters coupled to 
propagating free photons in ID (see Fig.l). We assume 
a small separation between emitters for simplicity, thus 
we can take an instantaneous (vacuum mediated) inter- 
action between emitters. It is so called Markovian limit, 
where the causal propagation time of photons between 
two emitters has been neglected [42]. Here we generalize 
the recently developed Bethe-ansatz approach including 
open boundary conditions to a chain of multiple emitters 
coupled to photons. We are able to derive the exact single 
and two-photon scattering states, and the corresponding 
photon transmissions in our full system. It is possible to 
extend our generalized approach to calculate scattering 
states and transmission for three or more photons. We 
here give a detailed calculation for a minimal model of 
the chain, namely a chain of two emitters, and briefly 
discuss results for a chain of three emitters in the next 
section. We take the following Hamiltonian to describe 
the two emitters-photons system. The full Hamiltonian 
H = Hq+Hs+Hc in real space within the rotating- wave 



A. Single-photon dynamics 

First we derive a single-photon scattering state and 
the corresponding transmission coefficient. We consider 
that an incident photon with wave-vector k (and energy 
Ef. = k) is injected from the left of the emitters. In 
our approach the wave-function at the left and the right 
sides of the emitters at x < describes the full system 
before scattering of photons from the emitters, and the 
wave-function at x > characterizes the system after 
scattering. The single-photon incident state is \ipl n ), and 
the wave-function after scattering is |"0out)- 

l^) = i/ dxe ik *a\(x)\0,l,l), 
l^out) = J dx[(<Pl(x)a\(x) + <Pl(x)al(x))\Q, 1, 1) 
+ ^)( e i|0,2,l} + 4|0,l,2))]. (2.2) 

Here \n, Zi,m 2 ) denotes state of the full system where n 
number of photons in the waveguide, left emitter in l\ 
state and right emitter in m 2 state, el is an occupation 
amplitude of the jth emitter in the excited state. We find 
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FIG. 2. Single photon transmission T k through two emitters 
versus energy E k of an incident photon. In all plots, Sl 2 = 
0.35, Ti = 0.01 and (a) «i = 0.35, T 2 = 0.01, J = 0.005, (b) 
Qi = 0.35, T 2 = 0.01, J = 0.05, (c) Qj_ = 0.2, T 2 = 0.01, 
J = 0.05 (blue full curve), fii = 0.35, T 2 = 0.0225, J = 0.05 
(maroon dash curve) and (d) fi a = 0.2, T 2 = 0.0225, J = 0.05. 



different amplitudes in IV'out) by solving the following lin- 
ear equations obtained from the stationary single-photon 
Schrodinger equation. 



-id x <t>l{x) - E k <f>l(x) + V 1 e\5{x) = 0, 
-id x 4> 2 k (x) - E k <&{x) + V 2 e 2 S(x) = 0, 
(fir - E k )e\ + V 1( l>i(x)6(x) + Je 2 = 0, 
[Q 2 - E k )el+V 2 cj>l{x)5{x) 



Je\ = 0. 



(2.3) 
(2.4) 
(2.5) 
(2.6) 



From Eqs.2.3,2.4 we get discontinuity of photon 
wave-functions across zero, cp k (0+) = <p k (0—) — 
iViel ,and <j> 2 k (0+) = <f>l(0-) - W 2 e\. Using the ini- 
tial conditions 4>\{Q—) = 1 and 0|(O— ) = for an inci- 
dent photon from the left, and the regularization across 

zero, 4(0) = [$(°+) + 4(°-)]A wc find <^(°) = 
1 -iVie\/2, 4> 2 k (0) = -iV 2 e 2 k /2 . These results along with 
Eqs.2.5,2.6 give the amplitudes of excited state, 
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where Ti = V, 



J 2 

E k -n 2 +iT 2 /2 



Vi 



(2.7) 



Next we calculate the amplitudes of outgoing photon 
states at the both sides of the emitters. 



X + ^x/2 

-iV X V 2 J 



{X + iT 1 /2)(E k -n 2 + iT2/2) 



e lkx 9(x), 



(2.8) 



where the single-photon transmission and reflection am- 
plitudes are respectively t\ = (x+tri /2X^-n 2 +*r 2 /2) and 



lFl / 2 For a same value of the left and the 
r, and identical emitters 



Cl, the transmission coefficient T k = \t\\ 2 



k x+iFi/2 

right coupling Y\ = T 2 

n x = n 2 

becomes one for E k — SI and J = Y/2 (see Fig. 2(a)), 
and the corresponding reflection coefficient is zero. We 
call it a single peak resonance (SPR). We plot transmis- 
sion coefficient through two emitters in Fig. 2 for different 
parameter sets. It shows that the transmission curve al- 
ways has two peaks except for a particular set of param- 
eters shown in Fig. 2(a). Also, the single-photon trans- 
mission curve becomes asymmetric in shape only when 
both Vtx 7^ f2 2 an d I\ 7^ T 2 . The transmission of a single 
photon can be detected by analyzing the temporal cor- 
relations of photons at the exit of the waveguide using 
single photon detectors for optical frequencies and linear 
detectors for microwave frequencies [45]. 



B. Two-photon dynamics 

It is easier to evaluate two-photon scattering state for a 
single emitter coupled to photons by introducing an even- 
odd transformation of the photons at the left and the 
right sides of the emitter. This transformation simplifies 
the calculation of the scattering state as the transformed 
even mode of photons is only coupled to the emitter while 
the transformed odd mode of photons is decoupled from 
the emitter [15, 18, 19]. Similar transformation of free 
photon modes is not useful for this problem as different 
emitters are coupled to photon modes at the left and 
right sides of the chain. Instead we derive two-photon 
scattering state directly using photon modes at the left 
and right sides of the chain. This path will also help us 
later to derive scattering states for a chain with multiple 
emitters. The two-photon incoming state |V>hi) f° r both 
the photons being injected from the left is given by, 



dx 1 dx 2 (/>k(xi,X2)'^=al(x 1 )a[(x 2 )\0, 1, 1), 



(2.9) 

where 4>k(x x ,x 2 ) = (e* 111 ^ 212 + e lklX2+lk2Xl )/2nV2 
with the incident wave vector k = (ki,k 2 ). The total 
energy of two incident photons E^ = ki + k 2 . We write a 
general two-photon scattering state IV'out) using the op- 
erators of free photon modes and emitters. 



= / dx\dx 2 



{gxx{xi,x 2 )^a\{xi)a\{x 2 ) 



+e{(xi)5(x 2 )a{(x 1 )<j 1+ + e\(xx)5{x 2 )a\ (xi)er 2 + 
+e 12 S(x 1 )S(x 2 )(T 1+ (T 2+ } + {gi 2 {xi;x 2 )a\(xi)a\{x 2 ) 
+e\{xi)8{x 2 )a\{xi)ai+ + e\{x{)5{x 2 )a\{x\)a 2+ } 

+g 22 {x 1 ,x 2 )-j=a\{x 1 )a\{x 2 )\ [0, 1, 1), (2.10) 

where gn(x x ,x 2 ) = g n (x 2 ,x i ) and g 22 (xi,x 2 ) = 
$22(2:2, x\) for the Bose statistics of photons. The ampli- 
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tudesgn(xi,x 2 ), 522 (#1, #2)1 and g 12 (x 1 ;x 2 ) denote out- 
going two-photon wave-functions, in which cither both 
photons arc transmitted or reflected, or one photon is 
transmitted while the other is reflected. A simple thought 
experimental set-up has been proposed in Rcf.[f5] to 
measure the magnitude of these wave-functions. In this 
set-up, a beam splitter with a single-photon counter on 
each of the output arm of the beam splitter, is placed 
at the entrance and the exit of the ID photonic waveg- 
uide. The measurement of the wave-functions are car- 
ried out by varying the distance of the photo-detectors 



after injecting a weak classical beam with the average 
photon numbers per pulse being much lower than two, 
and the pulse repetition rate being much smaller than 
the spontaneous emission rate of the emitters. Here 
e\(x), e\ (x) {e\{x) , e^ix)) are the probability amplitudes 
of one photon in the left (right) side of the emitters while 
respectively the left or the right emitter is in the excited 
state. Also, ei 2 is the probability amplitude of both the 
emitters are in excited state. We determine all the am- 
plitudes in Eq.(2.10) using the two-photon Schrodingcr 
equation. Thus, we obtain the following set of linear 
coupled first order differential equations. 



i{d Xl 


+ d X2 


) - E k ^gii(xi,x 2 ) + ^=[e\(xi)5(x 2 ) + 5(xi)e\(x 2 )] 


= 0, 


(2.11) 


(- 


id x - 


E k + fii) e\ (x) + [511 (0, x) + gil (x, 0)] + Jej (x) 


= 0, 


(2.12) 


(- 


- id x 


-E k + n 2 )e?(a;) + ^^(rr; 0) + Je\(x) + Vie 12 5(x) 


= 0, 


(2.13) 






(fii + fi 2 - E k )e 12 + Vie?(0) + V 2 e&0) 


= 0, 


(2.14) 


i(d Xl 


+ d X2 


) - #k) 522 (x i,x 2 ) + ^[el(x 1 )S(x 2 ) + S(x 1 )e 2 ,(x 2 )\ 


= 0, 


(2.15) 


(- 


~ id x 


-E k + Sl^4(x) + V l9l2 (0; x) + Jej(x) + V 2 e 12 5(x) 


= 0, 


(2.16) 






( y -id x -E k + n 2 )ei(ar) + ^2522(0,2;) + Je\{x) 


= 0, 


(2.17) 


i{d Xl 


+ d x . 


,) - E k ) g 12 (xi; x 2 ) + V 2 e\(x{)8{x 2 ) + Vi6{x{)e\{x 2 ) 


= 0. 


(2.18) 



The initial conditions for the amplitudes are chosen such 
that gu{x\,x 2 ), g 22 (xi,x 2 ), and g\ 2 (x\;x 2 ) at the re- 
gion x\,x 2 < correspond to plane wave- function of the 
Eq.(2.9). Thus, for both the photons being injected from 
the left, ffii (si < 0,2:2 < 0) = (j>k(x 1 ,x 2 ), g 22 (x x < 
0,x 2 < 0) = 0, 512(2:1 < 0;x 2 < 0) = 0. We use the fol- 
lowing regularization of these amplitudes across Xi,x 2 = 
as 511(0,2;) = gu(x,0) = [511(0+, x) + 5u(0-, x)]/2, 
g 22 (0,x) = g 22 ( Xl 0) = [522(0+, x) + g 22 (0-,x)}/2, 
5i2(0; x) = [512(0+; x) + 512(0-; x)]/2 and 512(2;; 0) = 
[512(2:; 0+)+ 512(2;; 0-)]/2. 

From Eqs. (2. 11, 2. 15, 2. 18) we get the discontinuity 
relations of the two-photon amplitudes across zero, 
g jj (x 1 ,0+) = g jj (x 1 ,0-) - iVj^(x{) / gjj (0+,x 2 ) = 

9jj(0-,x 2 ) - iV,e^x 2 )/V2, for j = 1,2, 



gi 2 {xi;0+) = gi 2 (xi; 0-) - iV 2 e\ (xi), 5i2(0+;2;2) = 
512(0— \x 2 ) — iV\e\{x 2 ). We also find the discontinu- 
ity relations of the emitter-photon amplitudes from 
Eqs.(2.12,2. 13, 2. 16,2.17), e}(0+) = eJ(O-), e^(0+) = 
e\(0-) - »Vie 12) 4(0+) = 4(0-) - iV 2 e 12 , e 2 2 (0+) = 
el(0-). Next we determine ei 2 from Eq.(2.14) and the 
previous relations for ef(0) and e 2 (0). We find 

P ^ief(0-) + ^(0-) . . 

12 n 1 + n 2 -E k -i(r 1 + v 2 )/2- 1 ' ' 

Now we derive the amplitudes of an excited emitter ex- 
plicitly. We rewrite the coupled linear inhomogeneous 
first-order differential equations in Eqs. 2. 12, 2. 13 in a ma- 
trix notation; we find 



-J \f e\(x) \ ( V2V l9ll (x,0-) \_( \ (n , 

h + l Vi/2 ) { e\(x) ) { V 29l2 (x- 0-) ) { V iei2 ) d ^' 2iJ > 



The eigenvalues of the above 2 x 2 square matrix (call A) \J a 2 + 16J 2 and a. = 2(Q 2 — Qi)+i(Ti — T 2 ). The inverse 
are A± = E k -(n 1 + n 2 )/2+i(r 1 +r 2 )/4 : ±f3/4:whcrc f3 = of the 2 x 2 square matrix V[= (|A_), |A+»] formed with 
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the eigenvectors |A_) and |A+) of A in Eq.2.20 is 



Therefore, we can write the coupled equations in 
Eq.2.20 in a rotated frame using the transformed exci- 
tation amplitudes of the emitters, e\(x) = jfe\(x) + 

°^e\{x) , e{{x) = - fe\{x) + =f^e?(a;) ) and the 

inverse transformations are ef(x) = e\(x) + ef(x) and 

e iW = lA - §) e ~iW " IA + The trans- 

formed matrix equation, 



* ( m )=>( v i ) ( m ) - *-* ( t^v ) - ( v,l ) '<*>■ 



We find the following discontinuity relations across zero 
for the transformed amplitudes from Eq.2.21, e}(0+) = 



ei(O-) 



V iei2 , e?(0+) = ef(O-) 



-a+i3 



Fie 



12- 



For the incident photons coming from the left 
512(2;; 0—) = at any x; we use the incident condition of 
511(2; < 0,0—) and find e\(x) and e\{x) at x < from 
the Eqs.2.21 



e\{x < 0) 
e\(x < 0) 

<Tfe 



7T/3 



(^e^+e^e**), where 

v_i 

k - (fii + si 2 )/2 + t(ri + r 2 )/4 - /3/4 ' 



J f 



+ Cfe 2 e* fcia: ), where 
Vi 



fc - (Oi + n 2 )/2 + i(ri + r 2 )/4 + /?/4 ' 



We now get back the original amplitudes for emitters' 
excitation using the inverse transformations. 



f 



el(x < 0) = ^(4^ + e 2 k y k n- (2-22) 



Physically these amplitudes describe part of the wave- 
function in which one photon is at x < of the incident 
photon channel while the other photon is absorbed by 
an emitter to reach to an excited state. The expres- 
sions in Eq.2.22 clearly reflect it. Similarly we rewrite 
Eqs.(2.16,2.17) in a matrix notation, and then transform 
it by multiplying I^ 1 from the left. Thus we find 



e\{x) 
e\(x) 







\(x) 



tt-uf V l9l2 (0-;x) 
\V2V 2 g 22 (0-,x) 



Viei2 




S(x) 



(2.23) 



where e\{x) - ^ 2 w t -3^ 
-^e\{x). As, <fe 2 (O-,aO,0i2(O-;a: 



for x < 0, 



we find e\{x < 0) = and e 2 (x < 0) = 0. Therefore, 
e\(x < 0) = e 2 (x < 0) = 0, i.e., the amplitude of one of 
the emitters being excited while one photon at x < at 
the right side of the emitters is zero. It is physically ac- 
ceptable, the causality prevents an injected photon from 
the left to move at x < of the right side photon channel. 
We now determine ei 2 , gu(x,0+), gn(0+,x), g 22 (x,0+), 
522(0+, x), gi 2 (x; 0+), and gi 2 (0+;x) at x < from the 
amplitudes of an excited emitter. We find explicit ex- 
pression for the amplitude ei 2 from Eq.2.19 by inserting 
ef(O-) and 4(0-). 



ei 2 



2nV 2 - sii - si 2 + i(ri + r 2 ) 



(2.24) 



which shows that ei 2 is non-zero only when a single pho- 
ton transmission amplitude is non-zero at least for one 



incident photon energy. The amplitude of one photon be- 
ing reflected while the other photon is still at the incident 
side is given by 

511 (x 1 < 0, x 2 > 0) = gu(x 2 > 0, x-i < 0) 

{r\ e ik 2 x i+ ik i x 2 _|_ r l e ik 1 x 1 +ik 2 x 2 -j (2 25) 



2\/27r 



Similarly we find the two-photon amplitudes, g 22 (xi < 
0,x 2 > 0) = ,g 22 (x 2 > 0,z;i < 0) = 0, 312(2:1 > 0;x 2 < 
0) = 0, and the amplitude of only one photon being trans- 
mitted while the other being at the incident channel is 



9i2 (2:1 < 0;a; 2 > 0) 
1 

~ 2^ 



^k 1 X 1 +ik 2 X 2t l^ +e ife 1 X 2 +ife 2 X lt lj^ (2.26) 



Next we calculate e\(x), e\{x), e\(x) and e 2 (a;) at x > 
from Eqs.2.21, 2. 23 using gu (x 1 < 0,x 2 > 0), g 22 (a;i < 
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0,x 2 > 0) and gi 2 (xi < 0;x 2 > 0). We find 



from Eq.2.26 and 322(0-, x) = 0. We get 



e\(x > 0) = Cl e iX -* + ^(4^^ + rl 2 e kl e ik >% 



e Y {x > 0) = c 2 e 



7T/3 



(r 1 c, p* 11 4- r 1 c, P ik 2 x \ 



where Ci and c 2 are two constants which we find using 
discontinuity relations of e\(x) and e\{x) across x = 0. 



ci = ^((1 - r-fejefe, + (1 - r fc2 ) £fel ) - i^_Fiei2, 
c 2 = --^((1 - rfejcfej + (1 - r k2 )<; kl ) - i Fiei 2 . 



2/? 



The general forms of e}(x) and ef(x) for any value of x 



are 



if \ f P f 1 a \ 
3l(a;)=; (2J ( 2-27 )Cl 



i\-x 



1 

2^ 



(4^1^) + 4 2 



e 2 (x) - (de^^ + ^e^^W 



(2.27) 



The first terms in e}(x) and e 2 (x) are contributions 
from a two-photon bound state. Similarly we calculate 
e\{x),e 2 (x) for x > from Eq.2.23 using 312(0— ; x > 0) 



~e\{x > 0) = cie"-* + ^(tle^* + tl 2 e kl e ik >% 



e 2 (x > 0) = 5 2 e jA ^ - -1(4^^ + i^e***), 

(2.28) 

where the constants c\ and c 2 are determined using dis- 
continuity relations of e\(x) and e 2 (x) across zero. 

J 2i J 

ci = -^p(t kl £k 2 +t k2 e kl ) - -yV 2 e 12 , 

J 2iJ 

C2 = ^(tfc^fe +*fc 2 ftj + -^-Vaei2. (2.29) 

One can find e 2 (x), e 2 (x) at any a; from e 2 (x), e 2 (x). 



4W=( 2 §4-^.e«--^(i + ^. W „ 



+ ^(4^(^ + 4^^))' 

s 2 (x) = (cie iA " 3; +C2e a + 3: )^(a;) 



(2.30) 



We are now ready to evaluate 
gu{xi,x 2 ), 922(^1, x 2 ), gi 2 (xi;x 2 ) for all values of 
xi,x 2 . Here we introduce a central of motion coordinate 
of photons x c = (xi + x 2 )/2 and a relative motion 
coordinate x = (xi — x 2 ). The general form of the 
two-photon amplitudes at different sides of the emitters 
are 



gu(xi,x 2 ) 



(<(Z1)^UZ2) + <^ 2 (Z1)<4>2)) 



iV x p ,1 a 



2^2/2 2/3 



)cie 



Ekic _*(A_-B k /2)|x| 



6{xi)6(x 2 ) 



+ ^}tA + §)^ < *"e 1(A+ -* /a)l " l »(«i)»(* a ). 



ff22 (x 1) x 2 ) = ^ 7 =(^ 1 (x 1 )^ 2 (x 2 ) + ^ 2 (x 1 )^ 1 (x 2 )) - ^(5 ie ^ e ^^-^/ 2 )N+5 2 e^e^^^/ 2 )N)0(x 1 )0(x 2 ) 



Wl ,(.n: r 2 ) = _(^ i(a; 0^ 2 (x2) + 4 2 (x 1 )^ l (x2))-^i(^c ie iA -^e i ^- A 1 



1 

x e j(£;k - A+);ri )6'(x 2 - xi)0(xi) - iy 2 (cie iA - Xl e i ( Bk - A -) X2 + c 2 e iX + Xl e^ Ek - x +^ x ')e(x 1 - x 2 )6(x 2 ). (2.31) 



The second terms in gn(xi,x 2 ), g 22 (xi,x 2 ) and 
<7i2(xi;x 2 ) are contributions from a two-photon bound 
state which arises due to coherent photon-photon interac- 
tions mediated by the emitters. We refer them as bound 
state as these terms fall rapidly with increasing \x\ — x 2 | 
or |x|. Two photons can exchange energy and momen- 
tum between themselves after scattering from the emit- 
ters due to resonant interactions at the localized region 
of the emitters. Though the total energy and momen- 



tum of the two incident photons remain the same after 
scattering. This energy-momentum distribution of pho- 
tons is also the origin of background fluorescence which 
can be conceived as an inelastic scattering of one photon 
from a composite transient object formed by the emitter 
absorbing the other photon. The strength of background 
fluorescence, i.e., inelastic scattering is much enhanced 
for the present system in a confined geometry compared 
to free space in 3D [46]. Here we define a total energy 
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FIG. 3. Two reflected photons wave-function \gn(xi,X2)\ 2 for one emitter (first column) and two emitters (middle and last 
columns) at various incident photon-pair energy detuning 5E and energy difference A. We use fli — f?2 = and Fi = T2 — T 
in all the plots for two emitters. The parameters are, (I) Top row: SE = A = 0, E kl — E k2 — Q — 0.35, F — 0.01, (II) Middle 
row: SE = 0.04, A = 0, E kl = E k2 = 0.35, Q = 0.33, T = 0.01, and (III) Bottom row: 8E = 0, A = -0.1, E kl = 0.3, E k2 = 
0.4, Q = 0.35, r = 0.01. The coupling J = 0.005 for the middle column and J = 0.05 for the last column. 



detuning SE = E^ — (Q1 + Q 2 ) and a energy difference 
A = E kl —E k2 . We discuss thoroughly various properties 
of the outgoing wave-functions by tuning 5E and A, and 
compare these results with those for a single emitter in a 
waveguide. 

First we consider form of these wave- functions at single 
photon resonance of the identical emitters for a same 
value of the left and the right couplings, i.e., E kl — E k2 = 
S7i = SI2 and V\ = V2. Around these parameter sets, the 
single photon transmission curve has a single peak for 
J = r/2 (check Fig. 2(a)). The wave-functions at this 
special point behave as 

Siifo.sa) = -_L e ^ e -rM/2 cos (r| x |/2), 

V 27T 

922(x 1 ,x 2 ) = --L e ^(l- e - r M/ 2 cos(r|z|/2)), 

V 2tt 

9i2(xi,x 2 ) = i e uSk*. e -rM/3 sin (r x /2). (2.32) 

7T 

These functions for a direct coupled two-level emitter- 
photons system in a ID waveguide are g~n(xi, x 2 ) — 



-e iBkX "e- r ^ 2 /(V27r), g 2 2(xi,x 2 ) = -e lEkX -(l - 
e - r \*\/ 2 )/(V2n), gu(xi;x 2 ) = - e ^ k x ce -r\x\/2 No _ 
tice that the form of the two-photon wave-functions here 
for a single emitter at resonance is different from that 
of the side-coupled emitter studied by Shen and Fan 
[15]. It is because Shen and Fan consider opposite signs 
for the group velocity of the left and the right moving 
propagating photon modes in the side-coupled model. 
In fact one will get back the results of Shen and Fan 
(within an overall minus sign for 322(^1,^2)) by chang- 
ing X\ — > —X\, x 2 — > — x 2 in our results for g 22 {xi,x 2 ), 
and x 2 — > —x 2 in our results for gi 2 {x\,x 2 ). Compar- 
ing two-photon wave-functions of a single and two emit- 
ters, we find that for the two emitters there is an extra 
sinusoidal oscillation to the exponentially decaying two 
reflected photon wave-function of a single direct coupled 
emitter. This interesting feature shows up because of 
our instaneous coupling between two emitters which is a 
valid approximation for small separation between emit- 
ters. A similar feature also persists in the two trans- 
mitted, and one transmitted-one reflected photons wave- 
functions. Though, when |i| is small, \gn(xi,x 2 )\ 2 and 
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FIG. 4. Two transmitted photons wave-function |<722 , ^2) | 2 for one emitter (first column) and two emitters (middle and last 
columns) at various incident photon-pair energy detuning 8E and energy difference A. We use £7i = SI2 = O and ri = T2 — T 
in all the plots for two emitters. The parameters are, (I) Top row: 5E = A = 0, E kl ~ Eu 2 = f2 = 0.35, V — 0.01, (II) Middle 
row: 8E = 0.04, A = 0, E kl = E k . 2 = 0.35, Q = 0.33, T = 0.01, and (III) Bottom row: 5E = 0, A = -0.1, E kl = 0.3, E k2 = 
0.4, Q = 0.35, r = 0.01. The coupling J = 0.005 for the middle column and J = 0.05 for the last column. 



\g22(xi,x 2 )\ 2 have similar features for a single and two 
emitters in a ID waveguide. The slow sinusoidal oscilla- 
tions along with a decay of wave-function with increasing 
distance separation |x| should manifest in the coherence 
measurement, such measurement of the g^ (r) func- 
tion in each case [15, 33, 45, 47]. Interestingly it has been 
shown recently that a normalized second-order correla- 
tion of stationary pulse intensity of a strongly interact- 
ing photon gas contains a cosine oscillating part which 
has been indicated as the Friedel oscillations characteriz- 
ing a strongly interacting gas [48]. We plot \gu(xi, x 2 )\ 2 
for a single and two emitters at single photon resonance 
in Figs. 3(a), (b), and \g 22 (xi, x 2 )\ 2 for a single and two 
emitters at single photon resonance in Figs. 4(a), (b). The 
oscillation due to the cosine term for two emitters is ev- 
ident in Fig. 4(b), and the corresponding curve for a sin- 
gle emitter has no oscillation. Though the oscillation 
for two emitters in the two reflected photons curve is 
not quite clear. \gi2(xi, ) 1 2 is shown in Fig. 5(a), (b) 
for a single and two emitters at single photon resonance. 
The sinusoidal oscillations for two emitters make a sharp 
change in the pattern of \gi2(xi,X2)\ 2 for two emitters 



compared to a single emitter. While \gi 2 (xi,X2 = xi)\ 2 
has a peak of height l/2ir 2 for a single emitter, it is zero 
for two emitters. Thus, one reflected and one transmit- 
ted photons for two emitters are anti-bunched at SPR. 
We show \gu(xi,x 2 )\ 2 , [522(211, x 2 )\ 2 and \g\2{x\, x 2 )\ 2 
for two emitters at J = 0.05 in Fig. 3(c), Fig. 4(c) and 
Fig. 5(c) which show a completely different behavior com- 
pared to the J = 0.005 case for two emitters or a sin- 
gle emitter. There are two peaks in the single photon 
transmission curve at J = 0.05 compared to one peak at 
J = 0.005. 

We plot |<7ii(xi, X2)| 2 and [522(2:1, x 2 )\ 2 for anon-zero 
total energy detuning and a degenerate incident energy 
of photons in Figs.3(d),(e),(f) and Figs.4(d),(e),(f). The 
overall features of two reflected and two transmitted pho- 
tons are qualitatively similar for a single and two emitters 
at 6E and A = except that there are more oscil- 
lations for two emitters. But the value of \gn(xi,x 2 )\ 2 
at x = in Fig. 3(e) for two emitters is higher compared 
to a single emitter. It signals a higher photon-photon 
correlations for two emitters than for a single emitter 
at a detuned total energy. Physically two emitters able 
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FIG. 5. One reflected and one transmitted photons wave-function | <7i2 (a^i , x 2)\ 2 for one emitter (first column) and two emitters 
(middle and last columns) at various incident photon-pair energy detuning SE and energy difference A. We use fii = Q2 — fi 
and Ti = T2 — T in all plots for two emitters. The parameters are, (I) Top row: SE = A = 0, Ek 1 = Ek 2 = ft = 0.35, V — 0.01, 
(II) Middle row: SE = 0.04, A = 0, E kl = E k2 = 0.35, fi = 0.33, T = 0.01, and (III) Bottom row: SE = 0, A = -0.1, E kl = 
0.3, E k2 = 0.4, fi = 0.35, T = 0.01. The coupling J = 0.005 for the middle column and J = 0.05 for the last column. 



to reflect photons more than a single emitter away from 
single photon resonance. We show the role of an en- 
ergy difference of two incident photons on behavior of 
\gn{x 1 ,x 2 )\ 2 and \g22{xi, x 2 )\ 2 in Figs. 3(g), (h),(i) and 
Figs. 4(g), (h),(i) where SE is kept zero. We find the be- 
havior of \gu(xi, X2)\ 2 for a single and two emitters at 
J = 0.005 (at SPR) is similar, though it has a very dif- 
ferent behavior away from SPR for two emitters. While 
Iffii(sci) X2)\ 2 oscillates with x for any value of x for the 
first case, it decays to zero with increasing x for two 
emitters away from SPR for SE = and A ^ 0. There is 
almost no transmission of two photons for , E k2 7^ ft 
at J = 0.005. The strength of instaneous coupling J 
between emitters completely changes the pattern of two 
transmitted photons correlation in Figs. 5(h), (i). The 
peak of |<7n(xi, X2)\ 2 at x = for an increasing A with 
SE = reduces to zero from its maximum value both 
for a single and two emitters. Therefore, the correlations 
of two reflected photons changes from bunching to anti- 
bunching, i.e., the emitters can induce an effective spa- 
tial attraction or repulsion between two reflected photons 
[15]. Here we find that we can even alternate the na- 



ture of photon-photon interactions of two reflected pho- 
tons in the two-emitters model by changing the coupling 
strength between emitters. It can be achieved by tuning 
the separation between emitters [36, 41]. For all values of 
SE and A, \g22(xi,X2)\ 2 =0 at x = which indicates that 
two transmitted photons are always anti-bunched at the 
forward direction for two emitters in a ID waveguide as 
they are for a single emitter in ID. The anti-bunching of 
two transmitted photons in our direct coupled emitters 
model occurs because two photons can not be emitted 
simultaneously by the right (second) emitter, and the 
transmitted photons in the right side of the waveguide 
are solely due emission from the right emitter. While the 
nature of |<7ii(xi, a^)) 2 and \gi2{xi,X2)\ 2 both for a single 
and two emitters is determined by interference between 
the incident and the emitted photons, thus they show 
various nature at x — depending on the parameters. 

We have shown features of \g\2{x\] X2)\ 2 away from 
SPR in the middle and the last rows of Fig. 5. While 
\gi2{xi,X2)\ 2 is symmetric as a function of x at SE 7^ 
and A = for T\ — T 2 in case of a single emitter, it is 
asymmetric as a function of x for two emitters at SE 7^ 
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and A = 0. The Fig. 5(e) shows that the overall feature of 
\gi2(x\; X2) 1 2 including a dip and a large peak is around 
X\ > X2- It indicates that reflected photon leaves the left 
emitter earlier than the transmitted photon leaves the 
right emitter. It is acceptable as the transmitted photon 
has to pass two emitters in the two-emitters model. We 
find 1512(^1; ) 1 2 is asymmetric as a function of x for 
both a single and two emitters at SPR when 5E = and 
A 7^ 0. In both cases a peak appears at x\ > x^ which 
again indicates that the reflected photon leaves the emit- 
ters before the transmitted photon at these parameter 
sets. We also find that the oscillations of \gyi(x\; a; 2 )| 2 
at large x for a single emitter disappear for two emitters 
when J = 0.05. Finally we show how stronger correlation 




X.|"X2 ^1 _ ^2 

FIG. 6. Two reflected photons wave-function |pxi(aci, ara)] a 
for two emitters at 5E — —0.15, A = 0, E kl — E k2 = Q 1 = 

0.35, fi 2 = 0.5, ri = r 2 = 0.01. 

between photons can be generated by cascaded optical 
nonlinearity of emitters when f2x 7^ 2 . The correlation 
between the two reflected photons for a degenerate inci- 
dent energy, i.e., A = at J = 0.005 and J = 0.05 are 
plotted in Fig. 6. We notice \gu(xi,X2 = xi)\ 2 = 0.454 at 
J = 0.005 is much larger than the maximum correlation 
\g%x(xi,X2 = xi)\ 2 = l/(27r 2 ) for a single emitter in a ID 
waveguide. 



C. Nonreciprocal photon transmission 

We define a photon current operator, I — i[H,N± — 
N 2 ]/2 where N± (N 2 ) is the total number of photons op- 
erator in the left (right) side of the emitter chain. There- 
fore, 7 = i(Vic\(G)a- - V 2 c\{Q)a- - H.c.)/2. We take 
expectation of current operator in the single and the two- 
photon scattering eigenstates of % to derive a single and 
two-photon current. The single-photon current for a pho- 
ton coming from the left (channel 1) is given by 

7(l;fc)H^tWoutH^l4| 2 - (2-33) 

Here we mention that 7(1; k) — 7(2; k) irrespective of 
the values of couplings and transition energies, including 
V\ 7^ V2 and fii 7^ 2 . Next we calculate photon cur- 
rent in the two-photon scattering states IV'out) f° r both 
the photons with wave- vectors k\ , k 2 being incident from 
the left (channel 1). The two-photon current 7(1, fci,^) 
has two parts, one 7o(fci,/c 2 ) is a contribution of two 



non-interacting photons, and the other <J7(1, fci, fo) is 
a change of two-photon current due to photon-photon 
correlation generated by resonant interactions between 
emitters and photons. The expression for 57(1, k\, k^) 
for general values of Q,j, Tj is quite long, and we have 
kept it in Appendix A. 

7(1, ki,k 2 ) — 7 (fci, k 2 ) + <57(1, ki,k 2 ), where 

7 (fci, k 2 ) = ^(\tl | 2 + |4 2 1 2 + 2|4 \ 2 5 klM )(2M) 

Here £ denotes length of the ID waveguide. The con- 
tribution in two-photon current from the non-interacting 
photons increases by a factor of two for two degenerate 
incident photons. This factor was left out in our pre- 
vious works [18-20]. We also notice that 51(1, k\, k 2 ) is 
one order of magnitude smaller than Io(k±, /C2). The non- 
interacting part of two-photon current is the same for 
forward and backward current for any value of emitter- 
photon coupling and transition energy of two emitters. 
The magnitude of 61(1, k\,k 2 ) is different for both the 
photons being incident from the left or the right side of 
the emitters when either Vi 7^ V2 or Oi 7^ O2 ° r V\ 7^ V2 
and Ox 7^ ^2- Instead of calculating forward 7(1, fci, A2) 
and backward 7(2, ki,k2) current separately to check 
asymmetry of two-photon current, we here demonstrate 
nonreciprocal optical transmission by interchanging the 
values of transition energy of two emitters in a calcula- 
tion for two-photon forward current. For example, we 
find 61(1, h,k 2 ) = -0.0001 for Q x = 0.5, Q 2 = 0.3, and 
<57(l,fci,fc 2 ) = -0.10017 for n± = 0.3, 2 = 0.5 while 
k x = k 2 = 0.3, ri = T 2 = 0.01 for both the cases. Thus 
we find that the total two-photon current as well as the 
interaction induced two-photon current change are differ- 
ent for an interchange of transition energy of two emitters 
keeping all other parameters the same, which is equiv- 
alent to say, |7(l,/ci,fc 2 )| 7^ |7(2, fc x , k 2 )\ for Oi 7^ 2 
even when Vi =V 2 . Previously we propose a few-photon 
optical diode [18] for a single emitter where it is essen- 
tial to have different left and right side emitter-photon 
coupling to get a nonreciprocal optical transmission via 
broken spatial inversion symmetry. The desired asym- 
metry in the emitter-photon coupling for a single emitter 
fixed at a position for a strong emitter-photon coupling 
is relatively difficult to achieve than to fix two different 
emitters in a ID waveguide. Thus our present proposal 
is much easier to implement in experiment. 

The mechanism behind diode like behavior in the two 
emitters is a redistribution of energy and momentum of 
photons after passing through the two emitters. The 
two-photon bound-state part of the wave-function is re- 
sponsible for the redistribution. The asymmetry in the 
transition energies generates an asymmetry in the dis- 
tribution of energy and momentum (via the two-photon 
bound state) of the photons coming from the left or the 
right of the emitters. For example, when two degenerate 
photons are injected at an incident energy on resonance 
with the first emitter, the photons can pass through the 
emitter, and their energy gets redistributed which helps 
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photons to pass the second emitter which has a different 
transition energy than the first. But if the first emit- 
ter's transition energy is different from the energy of the 
incident photons, photons can not transmit through the 
emitter and a redistribution of their energy does not oc- 
cur. Thus photons can not pass the second emitter in 
the last case even when its transition energy is the same 
as incident photon energy. Therefore, we get different 
contributions in the two-photon current change for the 
two cases. The multi-photon transport is also correlated 
in this system and the proposed system also acts as an 
optical diode for many photons. There are several dif- 
ferent proposal for realizing optical diodes using various 
mechanisms, such as, magneto-optic effect, macroscopic 
optical nonlinearity, indirect inter-band photonic transi- 
tions, opto-acoustic effect [49, 50]. Our proposed optical 
diode works at low intensity of light in the fully-quantum 
regime compared to most previous proposals in the classi- 
cal regime, and may have applications to build quantum 
logic gates for optical quantum information processing 
and quantum computation. The practical realization of 
the optical diode relies on achieving a strong emitter- 
photon coupling in the waveguide and a difference in the 
transition energy of the emitters. 



III. THREE EMITTERS IN A WAVEGUIDE 

In this section we briefly mention the main results of 
a single and two-photon scattering states for a chain of 
three identical emitters coupled to free propagating pho- 
tons in a ID waveguide. The full Hamiltonian of the 
emitters-photons system within the rotating-wave ap- 
proximation is H = Ho + H4 + He where, 

H = — i J dx [a\{x)d x ai{x) + a\{x)d x a 2 {x)], 

3 

H A = n|2)ij<2| + J(ai-a 2 + + <t 2 -(t 3+ + H.c), 

3=1 

H c - (Vcj 1+ai (0) + V<r 3+ a 2 (0) + H.c). (3.1) 

The part of the Hamiltonian H , H^ and Hp denote 
respectively the propagating photon modes, the emitters 
and the emitter-photon coupling. Here f2 is a transition 
energy of the identical emitters and J is a coupling be- 
tween emitters. V is a coupling strength of photons with 
the left and right emitters, a} (a;) and a 2 {x) are creation 
operators of a photon at the left and the right side of 
the emitters respectively. <jj- — c]gCj e (< 7 j+ = c j e c jg) * s 
again a lowering (raising) ladder operator of the jth emit- 
ter where Cj g {Cj e ) is a creation operator of the ground 
state |l)j (excited state \2)j) of the jth two-level emitter. 
In this section, we use bold character for those symbols 
which also appear before for the two emitters model. 



A. Single-photon dynamics 

The incident state \ipl D ) of a photon coming from the 
left of the emitter, and the single-photon wave-function 
after scattering (V'out) are 



l^out 



1 



2n 
1 



dx e ikx a\{x)\0, 1,1,1), 

dx[(<f>l(x)al(x) + 4>l(x)4(x))\Q, 1, 1, 1) 
+ S(x)(el\0, 2, 1, 1) + e||0, 1, 2, 1) + e||0, 1, 1, 2))] . 



The basis state \n,li,m 2 ,p3) denotes total n number of 
photons in the waveguide, the left most emitter (no.l) in 
l\ state, the middle emitter (no. 2) in m 2 state and the 
right most emitter (no. 3) in p% state. e 3 k is an amplitude 
of the jth excited emitter. The amplitudes of l^out) are 
derived by solving five linear coupled (differential) equa- 
tions obtained from the single-photon Schrodinger equa- 
tion. The method is similar to the one in Subsec.IIA. 
The amplitudes for an emitter being in the excited state 
are, 

i = V((E k -n + iT/2)(E k -n)- J 2 ) 

Gfc (E k - n + ir/2) 2 (E k -n) - 2J 2 (E k -n + ir/2) ' 

e 2 = i 

k {E k -Q + iT/2){E k -Q) -2.P 1 

3 _ J^V_ 

Gfc ~ {E k - Q + iT/2) 2 {E k -fl)- 2,P{E k -n + iT/2) ' 

(3.2) 

where T = V 2 , and the amplitudes of outgoing photon 
state at the both sides of the waveguide are 

4>\{x)=e lhx 6{-x)+rle ik * 6{x), 
d>l{x)=t\e ik *0{x), 



(3.3) 



where the single-photon transmission and reflection am- 
plitudes are respectively = — iVe\ and t\ = 1 —iVe},. 
The transmission coefficient becomes one at E k — f2, and 
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FIG. 7. Single photon transmission T,J through three identical 
emitters versus energy Ek of an incident photon. In both plots 
n = 0.35, r = 0.01 and (a) J = 0.0025, (b) J = 0.05. 

the corresponding reflection coefficient vanishes. We plot 
the single-photon transmission coefficient through 
three identical emitters in Fig. 7 for different values of 
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intcr-cmittcr coupling J. Single peak resonance arises 
near E k = £1 for a very weak coupling between emitters, 
i.e., J << T (see Fig. 7(a)) while there are three distinct 
resonant levels for a relatively stronger coupling between 
emitters. 



B. Two-photon dynamics 

We again consider that both the photons are injected 
from the left side of the emitters. The two-photon in- 
coming state is given by, 

iV'hi) = J dx 1 dx 2 (f> ii (x 1 ,x 2 )^a\(x 1 )a\(x 2 )\0, 1, 1(3)4) 



where k {% i,x 2 ) 



(e 



ikiXi J r ik2X2 



ik\X2+ik2X\ 



)/27iV2 



with k = (k\,k 2 ) and total energy = k± + k 2 . The 
outgoing two-photon scattering state is 



IV'out) 



{gii(xi,x 2 )-j=a\(xi)a\(x 2 ) + (e}(xi)oi + + e((xi)a 2+ + e 3 1 (x 1 )a 3+ )5(x 2 )a\(x 1 ) 



J dxidx 2 

+S(x 1 )S(x 2 )(e 12 a 1+ a 2+ + e 13 <r 1+ <7 3+ + e 23 a 2+ a 3+ } + {gi 2 (a;i; x 2 )a\(xi)a\(x 2 ) + (e 2 (a;i)<7i + + el(xi)(j 2+ 

+el(x 1 )a 3+ )S(x 2 )al(x 1 )} + g 22 (a;i, x 2 ) ^—a\(x\)a\(x 2 ) |0, 1,1,1). (3.5) 

V2 J 



e\(x) = 6(x)(- Cl e iXlX + c 2 e l ^ x + c 3 e lX3X ) 



The character of different amplitudes in IV'out) i s evident reflected) is at the left side of the emitters, 
from their definition. We obtain explicit expressions for 
these amplitudes by solving twelve linear coupled first- 
order differential equations derived from the two-photon 
Schrodinger equation. The procedure is similar to that 
in Subsec.IIB, but now we deal with two identical 3x3 
square matrices to find the emitter excitation amplitudes 
instead of 2 x 2 square matrices in Subsec.IIB. The am- 
plitudes for any two emitters in the excited state are given 

by 



+ ^ r (e 1 kl <t> 1 k2 (x)+e 1 k2 <t> 1 ki (x)), 



(3.6) 



a 
47 



4 J 4J' 

1 (e^LW+efe^W)- 



iT 
4J 
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Next we write the amplitudes of a single excited emitter 
when one photon (which is either at the incident side or 



where a = V32J 2 -T 2 , Ai = £ k -Q+f , A 2 = £ k -fi + 
ir/4-a/4, A 3 = E k -fl+iT/4+a/4:. The constants cj in 
the above relations are arising from the discontinuity of 
the amplitudes at x — due to emitter-photon resonant 
interactions, and they are, 
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where e k = V/(k — Cl + iT/2), <; k = V/(k — Q + iT/A — of a single excited emitter when there is one transmitted 
a/4), cp k = V/(k — Q, + iT /4 + a/4). The amplitudes photon at the right side of the emitters are 
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All the terms with Cj and Cj in the expressions of the 
excited emitters in Eqs.3.6-3.11 are contributions from 
the two-photon bound state which is the source of back- 
ground fluorescence in the system. Again we define a 



I 

central of motion coordinate x c = (x\ + x 2 )/2 and a rel- 
ative motion coordinate x = (xi — x 2 ) of photons. The 
amplitudes of two photons at different sides of the emit- 
ters are 



gll(*l,*2) = ^(<(si)0U*2) + ^(^ 



+ c 3 e 
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+ c 3 e tX3X2 e i{Ek - X3)xi )9(x 2 - x 1 )9{x l ) - iV {c x e lXlXl e i{E ^ Xl)x2 + c 2 e lX2Xl e^""^)* 2 
+ cse^V^-* 3 ^ 2 )^! - x 2 )6i(x 2 ). 



(3.12) 



(3.13) 



(3.14) 



The first parts of the expressions in Eqs.3. 12,3. 13,3. 14 
correspond to two non-interacting photons, while the 
rest of the terms involving the coefficients Cj and £j 
are a signature of strong photon-photon correlations. 
These correlations also represent background fluores- 
cence generated by inelastic scattering of a photon. We 
show these wave-functions in Figs. 8, 9 for different val- 
ues of a total energy detuning 8E = Ey — 2£1 and 
an energy difference A = Ej Cl — Ej. 2 . First, we plot 
|gn (x\,x 2 ) | 2 , |g 2 2(a;i,a;2)| 2 and |gi 2 (xi; x 2 )\ 2 at single 
photon resonance, Ek 1 = Ek 2 = SI in the first col- 
umn of Fig.8 for SE = A = and J = 0.005. While 
|gn(a;i,a;2)| 2 , ^22(^1, x 2)\ 2 are symmetric as a function 
of x, |gi 2 (xi; x 2 )\ 2 is always asymmetric including the 
single photon resonance. This is different from the single 
and the two emitters models. The peak in |gi 2 (xi; X2)| 2 
for J = 0.005 occurs at X\ > x 2 which indicates that 
the reflected photon leaves the first emitter before the 
emission of the transmitted photon from the right most 
emitter. The two transmitted photons are always anti- 



bunched for the three emitters for a same reason as stated 
for the two emitters. Notice that |gn(xi,x 2 = xi)| 2 = 
1/27T 2 at single photon resonance; this is similar to the 
single and the two emitters case. But now the form and 
magnitude of |gn(xi,x 2 = xi)| 2 , |g 2 2(zi, x 2 )| 2 for the 
three emitters are quite different from the single and 
the two emitters at single photon resonance. In fact we 
find a large magnitude for the peaks of |gn(a;i, £ 2 )| 2 , 
\g 22 (xi,x 2 )\ 2 and |gi2 (^1 ; ^2) | 2 at non-zero x. Next we 
show these wave-functions for a detuning of either SE or 
A in the second and third columns of Fig.8. We exam- 
ine that any detuning reduces the height of the peaks of 
these wave-function by large amount. Finally we show 
these wave-functions for a stronger emitter-emitter cou- 
pling in Fig. 9 for J — 0.05. An increase in the strength 
of coupling between the emitters increases the number 
oscillations in these wave- functions, and simultaneously 
reduces the height of the peaks. It can be understood 
from the single-photon transmission curve of the three 
emitters in Fig. 7 which shows that an increase in the 
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FIG. 8. Two transmitted photons wave-function |gn(xi, 2:2) | 2 , |<?22(2:i, £2)| 2 and |gi2(a;i, £2)| 2 for three identical emitters at 
various incident photon-pair energy detuning SE and energy difference A. In all plots J = 0.005. The parameters are, (I) First 
column: SE = A = 0, E kl = E k2 = SI = 0.35, T = 0.01, (II) Middle column: 5E = 0.04, A = 0, E kl = E k2 = 0.35, SI = 
0.33, T = 0.01, and (III) Last column: 5E = 0, A = -0.1, E kl = 0.3, E k2 = 0.4, tt = 0.35, T = 0.01. 



coupling increases the number of peaks in the single pho- 
ton transmission curve due to the separation of resonant 
energy levels. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we calculate an exact solution of the sin- 
gle and two-photon scattering states for free-propagating 
photons interacting resonantly with two and three two- 
level emitters in a ID photonic waveguide. The two- 
photon transport in these systems are strongly corre- 
lated because the two-level emitter can absorb only a 
single photon at a given time. We derive the single and 
two-photon transmission in the waveguide for two emit- 
ters, and show that the two-photon transmission is non- 
reciprocal when the two emitters have different transi- 
tion energies. Recently a large optical non-reciprocity 
is observed within telecommunication wavelengths in an 
experiment based on strong optical nonlinearity in two 
high-quality factor silicon micro- rings [50]. In our pro- 
posal we replace those micro-rings by two-level emitters, 



and employ a strong optical nonlinearity due to emitter- 
photon interaction in a confined geometry along with an 
inherent breaking of spatial inversion symmetry by the 
emitters with different transition energies. We expect 
that the proposed optical-diode for two different emit- 
ters can be realized easily using present experimental 
techniques. We have also studied how photon-photon 
correlations scale with an increasing number of emitters 
in the waveguide. Thus, our study provides a better un- 
derstanding to control coherent optical nonlinearity in 
small nanostructures. 

Here we generalize a recently developed fully quantum- 
mechanical transport approach for a single emitter to 
multiple emitters. The main feature of the approach is 
a proper inclusion of a two-photon bound state which 
represents background fluorescence due to inelastic scat- 
tering. The present generalized approach should be use- 
ful for a wide class of quantum transport problems in 
variety of systems including the mixed Bose-Fermi cold 
gases. Now we briefly list some areas where the method 
can be useful or further developed. First, the above ap- 
proach should be generalized further for three or more 
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FIG. 9. Two transmitted photons wave-function |3ii(xi,Xa)| 2 i 1322(2:1, X2)\ 2 and IP12 , ) | 2 for three identical emitters at 
various incident photon-pair energy detuning 8E and energy difference A. In all plots J = 0.05. The parameters are, (I) First 
column: SE = A = 0, E kl = E k2 = SI = 0.35, F = 0.01, (II) Middle column: 5E = 0.04, A = 0, E kl = E k ., = 0.35, SI = 
0.33, T = 0.01, and (III) Last column: 5E = 0, A = -0.1, E kl = 0.3, E k2 = 0.4, SI = 0.35, T = 0.01. 



photons. It has been already shown that the open-Bethe 
ansatz procedure works for multi-particle quantum trans- 
port with a single impurity [51]. It would be possible 
to go beyond single impurity for multi-photon case fol- 
lowing the generalized method of this paper. Second, 
it is also interesting to study the photon-photon cor- 
relations mediated by resonant interactions of photons 
with three or multi-level emitters. Recently, a scaling of 
electromagnetically induced transparency for weak light 
field with an increasing number of three-level emitters 
has been demonstrated in a beautiful experiment [10]. 
We can employ the technique of the present paper to 
study these on-going quantum optics experiments [52]. 
Third, it is equally relevant to have a general routine 
scheme of correlated two-photon transport for an arbi- 
trary number of emitters in ID. Fourth, one should also 
consider two or more channels of emitters in a quasi-lD 
geometry, and the resonant interactions of emitters be- 
tween intra and inter channels. The first step in this 
direction is to evaluate scattering states for two parallel 
emitters which are both coupled to the propagating pho- 



ton modes at the both sides of the emitters [52]. Fifth, 
numerical studies, such as using DMRG [22] or lattice 
Green's function [53] for discrete photonic lattices with 
finite band-edges should be extended for multi emitters, 
and a comparison of these numerical results with the 
exact results in continuum should be carried out [20]. 
Finally, progress in the above directions will enrich us 
enough to study non-equilibrium quantum phase transi- 
tion of light in the multi-particle emitter-photon systems 
in comparison with the proposed equilibrium quantum- 
phase transition of light [54] . 
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Appendix A: Photon-photon interaction induced two-photon current change for two emitters 
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